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The general construction of quasi-classically concentrated solutions to the Hartree type equation, based 



o 



o 
o 



Abstract 



on the complex WKB-Maslov method, is presented. The formal solutions of the Cauchy problem for this 



equation, asymptotic in small parameter h (h — » 0), are constructed with a power accuracy of 0(h N/2 ), where 
N is any natural number. In constructing the quasi-classically concentrated solutions, a set of Hamilton- 
Ehrenfest equations (equations for centered moments) is essentially used. The nonlinear superposition prin- 
ciple has been formulated for the class of quasi-classically concentrated solutions of Hartree type equations. 
i-C ■ The results obtained are exemplified by a one-dimensional Hartree type equation with a Gaussian potential. 
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Introduction 

The nonlinear Schrodinger equation 

{-id t + H(t,\V\ 2 )}V = 0, (0.1) 

where l~t(t, | V I / | 2 ) is a nonlinear operator, arises in describing a broad spectrum of physical phenomena. In 
statistical physics and quantum field theory, the generalized model of the evolution of bosons is described in terms 
of the second quantization formalism by the Schrodinger equation Q which, in Hartree's approximation, leads 
to the classical multidimensional Schrodinger equation with a nonlocal nonlinearity for one-particle functions, 
i.e., a Hartree type equation. 

The quantum effects associated with the propagation of an optical pulse in a nonlinear medium are also 
described in the second quantization formalism by the one-dimensional Schrodinger equation with a delta-shaped 
interaction potential. In this case, the Hartree approximation results in the classical nonlinear Schrodinger 
equation |2| |j| which is integrated by the Inverse Scattering Transform (1ST) method and has soliton solutions 
|Q. Solitons are localized wave packets propagating without distortion and interacting elastically in mutual 
collisions. The soliton theory has found wide application in various fields of nonlinear physics |5|, ^, |7], 

Investigations of the statistical properties of optical fields have led to the concept of compressed states of a 
field in which quantum fluctuations are minimized and the highest possible accuracy of optical measurements 
is achieved. The important problem of the correspondence between the stressed states describing the quantum 
properties of a radiation and the optical solitons is analyzed in J^, |!| . 

The Hartree type equation is nonintegrable by the 1ST method. Nevertheless, approximate solutions showing 
some properties characteristic of solitons can be constructed. Solutions of this type are referred to as solitary 
waves or "quasi-solitons" to differentiate them from the solitons (in the strict sense) arising in 1ST integrable 
models. 
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An efficient method for constructing solutions of this type is offered by the technique of quasi-classical 
asymptotics. Thus, for nonlinear operators of the self-consistent field type, the theory of canonical operators 
with a real phase has been constructed for the Cauchy problem in |], [l0| and for spectral problems, including 
those with singular potentials, in jn], [D| (see also ju], [H| Solitonlike solutions of the Hartree type equation 
and some types of interaction potentials have also been constructed in jl7] . 

In this paper, localized solutions asymptotical in small parameter h (h — > 0) for the (nonlinear) Hartree 



type equation are constructed using the so-called WKB method or the Maslov complex germ theory |18| ?, 21 
The constructed solutions are a generalization of the well-known quantum mechanical coherent and compressed 
states for linear equations [^4| |25| for the case of nonlinear Hartree type equations with variable coefficients. We 
refer to the corresponding asymptotical solutions, like in the linear case [^lj, as quasi- classically concentrated 
solutions (or states). 

Most typical of solitary waves ( "quasi-solitons" ) is that they show some properties characteristic of particles. 
For the "quasi-solitons" being quasi-classically concentrated states of a Hartree type equation these properties 
are represented by a dynamic set of ordinary differential equations for the "quantum" means X(t, h) and 
P(t, h) of the operators of coordinates x and momenta p and for the centered higher-order moments. In the 
limit of H — > 0, the centroid of such a quasi-soliton moves in the phase space along the trajectory of this 
dynamic system: at each point in time, the quasi-classically concentrated state is efficiently concentrated in the 
neighborhood of the point X(t, 0) (in the x representation) and in the neighborhood of the point P(t, 0) (in the 
p representation). Note that a similar set of equations in quantum means has been obtained in [19, |2(J for the 
linear case (Schrodinger equation) and in for a more general case. It has been shown 22| 23] that these 
equations are Poisson equations with respect to the (degenerate) nonlinear Dirac bracket. Therefore, we call 
the equations in quantum means for the Hartree type equation, like in the linear case |2l|] , Hamilton-Ehrenfcst 
equations. The Hamiltonian character of these equations is the subject of a special study. Nevertheless, it 
should be noted that, as distinct from the linear case, the construction of the quasi-classically concentrated 
states for the Hartree type equation essentially uses the solutions of the Hamilton-Ehrenfest equations. 

The specificity of the Hartree type equation, where nonlinear terms are only under the integral sign, is that 
it shows some properties inherent in linear equations. In particular, it has been demonstrated that for the 
class of quasi-classically concentrated solutions of this type of equation (with a given accuracy h, h — > 0), the 
nonlinear superposition principle is valid. 

In terms of the approach under consideration, the formal asymptotical solutions of the Cauchy problem for 
this equation and the evolution operator have been constructed in the class of trajectory-concentrated functions, 
allowing any accuracy in small parameter h, h — > 0. 

It should be stressed that throughout the paper we deal with the construction of the formal asymptotical 
solutions to the Hartree type equation with the residual whose norm has a small estimate in parameter fi,, h — ^ 0. 
To substantiate these asymptotics for finite times t £ [0, T], T = const, is a special nontrivial mathematical 
problem. This problem is concerned with obtaining a priori estimates for the solution of nonlinear equation 
(0.1), which are uniform in parameter h e]0, 1], and is beyond the scope of the present work. Note that, in view 
of the heuristic considerations given in |ll|, it seems that the difference between the exact and the constructed 
formal asymptotical solution can be found with the use of the method developed in Jll], |H] ■ 

This paper is arranged as follows: The first section gives principal notions and definitions. In the second 
section, a class of trajectory-concentrated functions is specified and the simplest properties of these functions are 
considered. In the third section, Hamilton-Ehrenfest equations are constructed which describe the "particle- like" 
properties of the quasi-classically concentrated solutions of the Hartree type equation. In the fourth section, the 
Hartree type equation is linearized for the solutions of the Hamilton-Ehrenfest equations, and a set of associated 
linear equations which determine the asymptotical solution of the starting problem is obtained. In the fifth 
section, we construct, accurate to 0(ft 3 / 2 ), quasi-classical coherent solutions to the Hartree type equation. In the 
sixth section, the principal term of the quasi-classical asymptotic of this equation is obtained in a class of quasi- 
classically concentrated functions. The quasi-classically concentrated solutions to the Hartree type equation 
are constructed with an arbitrary accuracy in \fh in Sec. 7, while the kernel of the evolution operator (Green 
function) of the Hartree type equation is constructed in Sec. 8. Herein, the nonlinear superposition principle 
is substantiated for the class of quasi-classically concentrated solutions. In the ninth section, the Hartree 
type equation with a Gaussian potential is considered as an example. Appendix A presents the properties 
of the solutions to a set of equations in variations necessary to construct the asymptotical solutions and the 
approximate evolution operator to the Hartree type equation. 
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1 The Hartree type equation 

In this paper, by the Hartree type equation is meant the equation 

{-ihd t + H(t) + xV(t, *)}* = 0, * e L a (K£). 

Here, the operators 

H(t)=H{z,t), 
V(t,*)= [ dy^*{y,t)V{z,w,t)-*{y,t) 



(1.1) 

(1.2) 
(1.3) 



are functions of the noncommuting operators 



x,y e 



the function iff* is complex conjugate to "J, x is a real parameter, and ft is a "small parameter", h E [0, 1[. For 
the operators z and it), the following commutative relations are valid: 



[zk,Zj]_ = [w k ,Wj]_ = ihJ kj , 
[Bk, Wj]_ = 0, k.j = 1, 2n, 



(1.4) 
(1.5) 



where J = || J k j\\2nx2n is a unit symplectic matrix 



J 



-I 

1 



For the functions of noncommuting variables, we use the Weyl ordering 27 
for instance, for the operator H.: 



In this case, we can write, 



(1.6) 



M 2 " 



where H(z, t) — H(p, x, t) is the Weyl symbol of the operator H(t) and (., .) is the Euclidean scalar product of 
the vectors 

n 2n 
(p, x) — ^^PjXj, p,xG R™, (z, w) = z j w j, Z,W £ K 2 ". 

i=i i=i 
We here are interested in localized solutions of equation ( |1 . l| ) for each fixed h € [0, 1[ and t € K, belonging 
to the Schwartz space with respect to the variable x £ W 1 . For the operators 7i(i) and V(t, iff) to be at work 
in this space, it is sufficient that their Weyl symbols TL(z,t) and V(z,w,t) be smooth functions^ and grow, 
together with their derivatives, with \z\ — > oo and |u>| — > oo no more rapidly as the polynomial and uniformly in 
t 6K 1 . Therefore, we believe that the following conditions for the functions Tt(z) and V(z,w,t) are satisfied: 

Supposition 1 For any multi-indices a, [3, [i, and v there exist constant Cp(T) and C^(T), such that the 
inequalities 



a 9^n(z,t) 



dzP 



< CUT), 



a u d^+ v \V(z,w,t) 

z - 



dz^dw v 



<C^(T), 



z,w <G 



< t < T 



are fulfilled. 



In what follows we assume that for all the operators under consideration, A = A(z,t), their Weyl symbols satisfy Supposition 
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Here, a, (3, /z, v are multi-indices (a, f3,fj,,v€ Z+") denned as 



a = (a 1 ,a 2 , ■ ■ ■ , a 2n ), M = ot\ + a 2 + ■ . . + a 2n , 



d^V(z) 
dz a 



dz? i dz?...dz%r 



cxj = 0, 00, j = 1, 2n. 



We are c omin g now to the description of the class of functions for which we shall find asymptotical solutions 
to equation (1.1). 



2 The class of trajectory-concentrated functions 



Let us introduce a class of functions singularly depending on a small parameter ft, which is a generalization of 
the notion of a solitary wave. It appears that asymptotical solutions to equation (1.1) can be constructed based 
on functions of this class, which depend on the phase trajectory z — Z(t, ft), the real function S(t, ft) (analogous 
to the classical action at x = in the linear case), and the parameter ft. For ft — > the functions of this class 
are concentrated in the neighborhood of a point moving along a given phase curve z = Z(t, 0). Functions of this 
type are well known in quantum mechanics. In particular, among these are coherent and "compressed" states 
of quantum systems with a quadric Hamiltonian [25| ^9|, |3(], [n], g3|, |34|, [55| ||(| . Note that the soliton 
solution localized only with respect to spatial (but not momentum) variables does not belong to this class. 
Let us denote this class of functions as V\{Z(t, H), S(t, H)) and define it as 



vl = vl(z{t,h),s{t,hj) 

$ : $(x,t,ft) 



/Ax , 

^Wh exp 



-(S(t,h) + (P(t,h),Ax)) 



(2.1) 



where the function y>(£, t, K) belongs to the Schwartz space § in variable £ S K™ and depends smoothly on t and 
regularly on \fh for h — > 0. Here, Ax = x — X(t, and the real function S(t, H) and the 2n-dimensional vector 
function Z(t, h) — (P(t, h),X(t, ft)), which characterize the class Vl(Z(t, h),S(t, ft)), depend regularly on vft in 
the neighborhood of ft = and are to be determined. In the cases where this does not give rise to ambiguity, 
we use a shorthand symbol of for V^(Z(t, ft), S(t, ft)). 
The functions of the class are normalized to 



\mf = (m\m) 



in the space L 2 (M.™) with the scalar product 



(*(*)!*(*)) 



In the subsequent manipulation, the argument t in the expression for the norm may be omitted: || < I > (t)|| 2 = 
||$|| 2 . 

In constructing asymptotical solutions, it is useful to define, along with the class of functions V\{Z{t, ft), Sit, ft)), 
the following class of functions 



ci{z{t,h),s{t,h)) : 

$ : $(x,t,h) 



/Ax . 
^tjexp 



;(S(t, K) + (P(t, ft),Af)) 



(2.2) 



where the functions if, as distinct from (2.1), are independent of ft 



At any fixed point in time ( e I 1 , the functions belonging to the class are concentrated, in the limit 
of ft — > 0, in the neighborhood of a point lyi ng o n the phase curve z — Z(t,0), t E M. 1 (the sense of this 
property is established exactly in theorems 2^-2^ below). Therefore, it is natural to refer to the functions of 
the class P\ as trajectory- concentrated junctions. The definition of the class of trajectory-concentrated functions 
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includes the phase trajectory Z(t, K) and the scalar function S(t, H) as free "parameters" . It appears that these 
"parameters" are determined unambiguousl y fro m the Hamilton-Ehrenfest equations (see Sec. 3) fitting the 
nonlinear (xr ^ 0) Hamiltonian of equation (1.1). Note that for a linear Schrodinger equation, in the limiting 
case of h — 0, the principal term of the series in h . — > determines the phase trajectory of the Hamilton 
system with the Hamiltonian H(p,x,t), and the function S(t, 0) is the classical action along this trajectory. In 
particular, in this case, the class Vt includes the well-known dynamic (compressed) coherent states of quantum 
systems with quadric Hamiltonians when the amplitude of tp in (|]l]) is taken as a Gaussian exponential: 

V {t,t) = exp[^{lQ(t)0]f(t), 

where Q(t) is a complex symmetrical matrix with a positive imaginary part, and the time factor is given by 

t 

Im Q(t) dr 



f(t) = 0m g(i)exp 



(see for details pi[). 

Let us consider the principal properties of the functions of the class V\{Z{t, H), S(t, H)), which are also valid 
for those of the class C\{Z{t, K), S(t, h)). 

Theorem 2.1 For the functions of the class V^(Z(t,H), S(t,h)), the following asymptotical estimates are valid 

($|{Az} Q |$) 



for centered moments A a (t, K) of order \a\, a € Zj™ • 



A a (t,h) 



(2.3) 



Here, {Az} a denotes the operator with the Weyl symbol (Az) a , 

Az = z - Z(t,h) = (Ap, Ax), Ap = p- P(t,h), Ax=x- X(t,h). 
Proof. The operator symbol {Az} a can be written as 

(Az) a = (Af) a * {Ax) a * , (a p ,a x ) = a, 



and, hence, according to the definition of Weyl-ordered pseudodifferential operators (1.6), we have for the mean 
value a a (t, H) of the operator {Az} a : 



($\{Az} a \<S>) = J dxdydp<S>*(x,t,h)x 



x exp 



(^((x-y),p))[Ap\^ 



Ax + Ay 



Ay = y-X(t,h). 



a a (t,H) 

Here, we have denoted 
After the change of variables 

Ax = Vh£, Ay = VhC, Ap = Vhio 
and taking into consideration the implicit form of the functions 

$>(x,t,H) = exp{i/K(S(t,H) + (P(t,n),Ax))}ipf^,t,H 

belonging to the class V\(Z{t, h), S(t, h)), we find 

<r a (t,h) = -J-^h^h^-M [ d£dtdujy*{it,h)x 

(2-kK) J 
x exp{»<(£- C), uj)}u} a * (f + CTMC, t, h) = 

= a ( " +|a|)/2 M a (t, n), 
H$l|2 = a"/ 2 [d£p*(g,t,h)<p(g,t,h) = h n / 2 M (t,h). 



(2.4) 
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Since <p(£, t, h) depends on y/h regularly and Mo(t, h) > 0, we get 



Aa ^ h >- ||$||2 - h Mo (t,h) * 

<fi M/a max =0(^1/2), 
te[o,r] M (t,H) v 7 

and thus the theorem is proved. 

Let us denote by the symbol 0(h v ) an operator F, such that for any function $ belonging to the space 
V^(z(t,H),S(t,h)) the asymptotical estimate 



1^1 
11*11 



o(nn, 



ft-> o, 



is valid. 



Theorem 2.2 For the functions belonging to T ,t h (Z(t 1 H), S(t, H)), the following asymptotical estimates are valid: 



{Az} a = 6(fr H/2 ) 
Proof is similar to that of relation (pT3|). 



(2.5) 



Corollary 2.2.1 For the functions belonging to V^(Z(t,h), S(t,h)), the following asymptotical estimates are 
valid: 



{-ihd t - S(t,h) + (P(t,h),x(t,h)) + (z(t,h),JAz)} = 6(h), 

Ax k = 6(Vh), Apj = 6(Vh), k,j = T~n. 



(2.6) 
(2.7) 



Proof follows from the explicit form ( }2.4| ) of the trajectory-concentrated functions K) G (2.2)] and 

from the estimates (2.5). 



Theorem 2.3 For any function $(ir, t, h) E V^(Z(t, h), S(t, h)), the limiting relations 

1 !<!>(£,£, ft) | 2 = S(x-X(t,0)), 



lim ,, 

h^o ||$|| 2 

1 



lim i,S)| 2 = 5(p - P(t, 0)), 

ft— >0 $ M 

where <&(p,t,h) = F%^^^[x^t,K), F^g^p is the frr 1 Fourier transform pq| , are uaM. 

Proof. Let us consider an arbitrary function 4>(x) G S. Then for any function $(x, £, h) G T 3 ^ the integral 



(2.8) 
(2.9) 



|#(t,a)| s 



|*(t,n)|| 



l$(t,< 



(x)\$(x,t,h)\ 2 dx = 



\<P(tM\ 



/Ax 



dec 



after the change of variables £ = Ax/Vh, becomes 

h n/2 



(|*(t,ft)| 2 0> 



0(x(t,R)+v^Ki;* > fi)i 2 df. 



Let us pass in the last equality to the limit of h —> 0, and, in view of 

M^)|| 2 = ft"/ 2 / i^ft)! 2 ^ 



and a regular dependence of the function ip(£,t, K) on yh, we arrive at the required statement. 

The proof of relation (2.9) is similar to the previous one if we notice that the Fourier transform of the 
function $(:?, t, K) G V\ can be represented as 



where 



*(£ t, K) = exp{ l - [S(t, h) - (p, X(t, h))] } 

<p{ij,t,h) = 



p-P{t,h) 



(2tt)™/ 2 



,t,h 



Denote by (L(t)) the mean value of the operator L(t), t S R 1 , self-conjugate in L2(M^), calculated from the 
function t, K) G V\. Then the following corollary is valid: 

Corollary 2.3.2 For any function t, h) G V\(Z(t, ft), S(t, h)) and any operator A(t, ft) whose Weyl symbol 
A(z,t,h) satisfies Supposition the equality 

lwa(A(t,ft)) = )im 1 t^{$($,t,H)\A(t,H)\$(x,t,H)) = 



ft-»o 



r,~ o <E> 



A(Z(t,0),f,0) 



(2.10) 



is valid. 



Proof is similar to that of relations ( |2.8| ) and (2.S). 
Following we introduce 



Definition 2.1 We refer to the solution <&(:?, t, ft) G V\ of equation (l.T) as qua si- classically concentrated 



the phase trajectory Z(t, H) for H — > 0, provided that the conditions (2.1) and t\2.9j ) are fulfilled 



Remark 2.1 The above estimates (2.5) of operators \Az } a allow a consistent expansion of the functions of 
the class V\(Z{t, H), S(t, h)) and the operator of equation (1.1) in a power series for y/h. This expansion gives 
rise to a set of recurrent equations which determine the sought- for asymptotical solution of equation (1.1). 
For any function $ G Vt(Z(t, ft), S(t, ft)), the representation 



N 



$(x, t,h) = J2 ft fc/2 $ (fc) (x, t, ft) + o(ti N+1 v 2 ), 



(2.11) 



k=0 



where ^ k ^(x,t,h) G C^(Z{t,ft), S(t,ft)), is valid. Representation ( 2.11 ) naturally induces the expansion of the 
space V\(Z{t, K), Sit, ft)) in a direct sum of subspaces 



P t h (Z(t,h),S(t,ft)) = ® Vi{Z{t,K),S{t,h),l). 



1=0 



(2.12) 



Here, the functions $ G T > t(Z(t,h), S(t,ft),l) C T > t(Z(t,h),S(t,h)), according to (2.2), have estimates by the 
norm 



1 



h n/2 



(2.13) 



where the function /i(t) is independent o f h a nd co ntin uously differentiable with respect to t. 
Similar to the proof of the estimates ( |2.5| ) and (|2.6|) ; it can be shown that the operators 

{Az} a , {-ihdt - S(t, h) + (P(t, ft),X(t, h)) + (Z(t, ft), JAz)} 
do not disrupt the structure of the expansion ( |2.11 ), ( 2.12| ), and 

{Az} a : V l h (Z(t, ft), S(t, h),l) -> Vl(Z(t, h),S(t, h), I + \a\), 

{-ihdt - S(t, K) + (P(t, K), X(t, h)) + (Z(t, ft), JAz)} : 
: V\(Z{t, h),S(t, h), I) -» vi(z(t, ft), S(t, K),l + 2). 



(2.14) 
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Remark 2.2 From Corollary [2.3.2 it follows that the solution ^(x,t, ft) of equation (fTl|), belonging to the class 



V\„ is quasi- classically concentrated. 



The limiting character of the conditions (2J5) and (2.9) and the asymptotical character of the estimates 



(2.3)-(2.6) valid for the class of trajectory-concentrated functions make it possible to construct quasi-classically 
concentrated solutions to the Hartree type equation not exactly, but approximately. In this case, the Li norm 
of the error has an order of ft", a > 1 for ft — > on any finite time interval [0, T]. Denote such an approximate 
solution as ^ as = ^as(^ t, ft)- This solution satisfies the following problem: 



•ft 9 
- lH di 



■H(t)+xV(t, * a 



0(ft"), 



$ a9 6^(U),S(t,S),fi), te [0,T], 



(2.15) 
(2.16) 



where 0(h a ) denotes the function g( a >(x,t,H), the "residual' of equation (1.1). For the residual, the following 
estimate is valid: 

U a \x,t,h)\\=0(h a ), ft->0. (2.17) 



max 

0<t<T ' 



Below we refer to the function \& as (x, t, ft) satisfying the prob lem ( |2.15|) -( 2.17 ) as a quasi-classically concentrated 
solution (modft", ft — * 0) of the Hartree type equation (p~l|). 

T he m ain goal of this work is to construct quasi-classically concentrated solutions to the Hartree type equa- 
tion (Q with any deg ree o f a ccura cy in small parameter y/h, ft — > 0, i.e., functions 4 , a s(^, t, ft) = \P^(af, t, ft) 
satisfying the problem ( 2.15 )— ( 2.17 ) in mod(ft' JV+1 ' / ' 2 ), where N > 2 is any natural number. 

Thus, the quasi-classically concentrated solutions ty( N \x,t,h) of the Hartree type equation describe ap- 
proximately the evolution of the initial state ^0(2^1 ft) if the latter has been taken from a class of trajectory- 
concentrated functions V\. The operators Ti.(t) and V(t,i&) entering in the Hartree type equation (1.1) leave 
the class V\ invariant on a finite time interval < t < T since their symbols satisfy Supposition [lj. Therefore, in 
constructing quasi-classically concentrated solutions to the Cauchy problem, the initial conditions can be taken 
in the form 

*(2,t,h)\ f= o = * D (0,h), * eV%(z ,S ). (2.18) 
The functions from the class V 1 ^ have the following form: 

*„(2, &) = ex P {^[S(0, ft) + (Po(ft), (x - Mh)))}}vo( X ~^ ih) ,h) , <Po& h ) e S(R£), (2.19) 

where Zq(K) — (Po(ft),Xo(ft)) is an arbitrary point of the phase space Rp™, and the constant So (ft) can be put 
equal to zero without loss of generality. 

Important particular cases of the initial conditions of type ( [2.19] ) are 



1) 



MO 



- „-<eUD/2 



where the real nxn matrix A is positive definite and symmetrical. Then relationship ( [2.19] ) defines the Go 
packet; 



2) 



where the complex nxn matrix Q is symmetrical and has a positive definite imaginary part Im Q and H u (if) and 
if G R™ are multidimensional Hermite polynomials of multi-index v = (v\, . . . , v n ) Jl2| . In this case, relationship 
(2.19) defines the Fock states of a multid ime nsi onal oscillator. 

The solution of the Cauchy problem (|l . l|) , ( [2.18 ) leads in turn to a set of Hamilton-Ehrenfest equations to 
the study of which we are coming. 



3 The set of Hamilton-Ehrenfest equations 

In view of Supposition [l] for the symbols TL(z,t) and V(z,w,t), the operator H(z,t) ( |l.2| ) is self-conjugate to 
the scalar product (^|$) in the space and the operator V(z,w,t) ( |l~3| ) is self-conjugate to the scalar 



8 



product L 2 (M 2 X y) 



*(t)|$(t)) R 2» = / dxdy**(x,y,t,h)${x,y,t,h). 



Therefore, for the exact solutions of equation (LI) we have 



||*(i)|| 2 = ||*(0)|| 2 , 

and for the mean values of the operator A(t) = A(z, t), calculated for these solutions, the equality 

d /dA(t)\ i 



H dy^*(y,t,h)[A(t),V(z,w,t)}^(y,t,h)y (3.1) 



where [A, £>]_ = AB — BA is the commutator of the operators A and B, is valid. We refer to equation (3.1) 
as the Ehrenfest equation for the operator A and function ^(x, t, ft). This term was chosen in view of the fact 
that in the linear case (x = 0) equation fll.l|) goes into the quantum mechanical Schrodinger equation and 



relationship (3.1) into the Ehrenfest equation ]26| 



Let us make the following notations: 

z=(fj,~x), Z(t,K) = (P(t,h),X{t,K)), Az = z-Z(t,h). (3.2) 
Using the rules of composition for Weyl symbols |27[] , we find for the symbol of the operator C = A B 

cW-A(3+2jA)flW-fl(J-»4K)- («> 

Here, the index over an operator symbol specifies the turn of its action. We suppose that for the Hartree type 
equation (JO]), exact solutions [or solutions differing from exact ones by a quantity 0(H°°)} exist in the class of 



trajectory-concentrated functions. Let us write Ehrenfest equations (3.1) for the mean values of the operators 
Zj,{Az} a calculated from such (trajectory-coherent) solutions of equation (1.1). After cumbersome, but not 
complicated calculations similar to the calculations that were performed for the linear case with x = (see 
for details we then obtain, restricting ourselves to the moments of order N, the following set of ordinary 

differential equations: 

N 1 N 1 

z= — J ( w *M(*>*)A M + >r -V zliV {z,t)A l ?j, 

\m\=o ^' kl=o ' 
Aa = ,^ ^^ 7 !(a-7)!(0-7W ( } 

N ^ 2n 
X ^H^(z,t) + k ^ —y^( z ,t)^Sj^a--/ + Jfj-Jj - ^ Zk<Xk A a (k) 
\v\=0 ' fc=l 

with initial conditions 



*lt=o = z ° = <*o|*|*o>, A Q | t=Q = (* |{z-z } Q |*o), (3.5) 
a £ Z 2 ", |a| < N. 



Here, iir = x||^oG?j h)\\ and ^(x, ft) is the initial function from ( 2.1 



d^H(z,t) d^+ u \V{z,w,t) 



H Zj[i (t, ft) = d z Hv(t, ft), a=(ap,a a: ), Ja = (oa, a p ), (3.6) 

2n 

(a - /3) = n ^("fc - Afc)> a(fc) = (ai - . . . , a 2 „ - S 2n ,k)- 

k=l 
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By analogy with the linear theory (x = 0) |2l[| , we refer to equations (3.4) as Hamilton-Ehrenfest equations 
of order N. In view of the estimates (^^), these equations are equivalent, for the class 7-^, to the nonlinear 
Hartree type equation (1.1) accurate to 0(f^ N+1 ^ 2 ). 

For the case of N = 2, the Hamilton Ehrcnfcst equations take the form 



where 



z = Jd z (l + ~(d z ,A 2 d z 

+%V(z,U,t))\ u)=z , 

[ A 2 = JMA 2 - A 2 MJ, 



M = [H zz (z,t) + xV zz (z,uj,t)] 



-{d u ,A 2 d u ))(H{z,t)- 



Equations (3.7) can be written in the equivalent form if in the second equation we put 

A 2 (t) = A(t)A 2 (0)A+(t), 

and then it becomes 



A = J MA 



A(Q) = I. 



(3.7) 



(3.8) 



4 Linearization of the Hartree type equation 



Let us now construct a quasi-classically concentrated (for h — > 0) solution to equation (1.1), satisfying the initial 
condition (2.18). 
Designate by 



yWfah) = (%,A^. 3 ,A^. 5i6 ,...) = (Z(t,h),A a (t,h)), \a\ < N 



,(2) A (3) 



(4.1) 



the solution of the Hamilton-Ehrenfes t equ ations of order N (^J) with the initial data y^ (0, K) (|3.5| ) deter- 
mined by the initial function ^o(x, H) ( 2.18 ), i.e., the mean values Z(0, H) and A Q (0, K) are calculated from the 
function "to (a;, K). Let us expand the "kernel" of the operator V(t, \E") in a Taylor power series for the operators 
Aw = w - Z(t,H): 

V{z,w,t) = 



\a\=0 



ai 



dw a 



=z(t,h) 



{Aw} c 



(4.2) 



Substituting this series into equation (1.1), we obtain for the functions 'f S V l n 



-ihdt + H(z, t) 



. ^ 1 d^V(z,w,t) 



a 



q|=0 

0(^+1)72^ 



dw a 



t=0 



=z(t,h) 



A a (t, h) 



* = 



where 



Z(t, h) = 

A a (t,h) = 



1 



\*(t,h)\\ 2 
l 



|*(t,i 



(*(t,h)\z\*(t,h)) ; 
5 (*(t,fi)|{Ai}«|*(i,n)). 



(4.3) 



(4.4) 



In view of the asymptotical estimates ( |2.3| ), the functi ons z(t, h) and A a (t, fi) can be determined with any degree 
of accuracy from the Hamilton-Ehrenfest equations (|3.4| ) as 



z(t, K) = z(t, h, N) + 0(h( N +v/ 2 ): 

A a (t, K) = A a (t, h, N) + 0{ti N+1 )/ 2 ) \a\ < N, 



(4.5) 



where z(t,H,N) and A a (t,H, N) are solutions of the Hamilton-Ehrenfest equations of order N, which are 
completely determined by the initial condition of the Cauchy problem for the Hartree type equation, ^o(x, t, H), 
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and do not use the explicit form of the solution ^5>(x,t,h) (4.3), Thus, the c hang e of the mean values of the 
operators for the solutions of the Hamilton-Ehrenfest equations of order N (4.5) linearizes the Hartree type 
equation (4.3) accurate to 0(h^ N+1 ^ 2 ). So, to find an asymptotical solution to the Hartree type equation ([l.lj ), 
we should consider the linear Schrodinger type equation: 



U N \t, * )$ - 0(ti N +^/ 2 ), $| t=0 = 4> ; 



(4.6) 



L W (t, ) = -ifift + H% t) + %J2 



N 1 d^V(z,w,t) 



|ct|=0 



dw a 



w=Z(t,h,N) 



A a (t, ft, JV). 



Definition 4.2 We caZ/ an equation of type ( |4.6|) wi/i a given a Hartree equation in the trajector y- co herent 



approximation or a linear associated Schrodinger equation of order N for the Hartree type equation ( 1.1 

The following statement is valid: 

Statement 4.1 If the function <&( N ^(x, t, h, ^o) € ^ *s an asymptotical (accurate to 0(tv- N+i ^ 2 ), H — > 0) 
solution of equation (4.6), satisfying the initial condition $|t=o = the function 



an asymptotical (accurate to 0(h^ N+1 ^ 2 ), h — > 0) solution of the Hartree type equation ( p..lj ). 
Now we expand the operators 



W(z,t), 



<9HV(z,w,i) 



w=Z(t,h,N) 



in a Taylor power series for the operator Az and present the operator —ihd t in the form 

-ihd t = {-(P(t,h,N),X(t,h,N)) + S(t,h)} - (Z(t,h,N),JAz) + 
+{-ihd t - S(t, h) + (P(t, h, AT), X(t, h, N)) + (Z(t, h, N), JAz)}. 



Here, the group of terms in braces containing —ihdt, in view of (2.14), has an order of 0(K). Other terms can be 
estimated, in view of ( |2.6| ), by the parameter h. Substitute the obtained expansions into ( |4.6| ). Take (accurate 
to 0(h N l 2 )) the real function S(t, K) entering in the definition of the class Vt(Z(t, h), S(t, K)) in the form 



S(t,H) = S^(t,H) = / \(P(t,h,N)X(t,h,N)) -H(Z(t,h,N),t)- 



o 

A 1 d^V(Z(t,h,N),w,t) 



|ee|=0 



a 



dw a 



w=Z(t,h,N) 



dt. 



(4.7) 



As a result, equation (4.6) will not contain operators of multiplication by functions depending only on t and h. 



In view of the estimates (EJJ) and (2.6) valid for the class V^(Z(t, h), S(t, h)), we obtain for (4.3) 



with the following notations: 



N 



h {N) (t^o) = Y, hk/2 ^( t ^o), 



fc=l 



Sj (t, * ) - S(t, h) + (P(t, h), X(t, h)) + (Z(t, h), JAz) + -{Az,Sj zz (t, * )A5), 
fizz (t, *o) = \n zz (z, t) + kV zz (z, w, t)] 

L J z=w=Z(t,h,N) 



(4.8) 

(4.9) 
(4.10) 
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' Ctl OZ r—'S.tl h M\ 



|a|=fe+2 



z=Z(t,h,N) 



E 



\a+f3\=k+2 



1 d\ a +P\y( z ,w,t) \ 
a!/3! dzPdw a 



\z=w=Z(t,h,N) 



{MYA a {t,h,N). 



(4.11) 



Here, k = 1,7V and the functions Z^(t) are the coefficients of the expansion of the projection Z(t,h) of the 
solution y( N >(t, H) of the Hamilton-Ehrenfest equations on the phase space R 2 ™ in a power series of the regular 
perturbation theory for \PK: 



N 



Z(t, h) = Z(t, h, N) = Z(t, o) + ft k/2 Z( k ) (t). 



k=2 



From the Hamilton-Ehrenfest equations, in view of the fact that the first-order moments are zero (A Q (i, K, N) = 
for |a| = 1), it follows that the coefficient Z^{t) is equal to zero. 

Remark 4.1 The solutions of the set of Hamilton-Ehrenfest equations depend on the index N that denotes the 
highest order of the centered moments A a , a € Z^™. We shall omit the index N if this does not give rise no 
ambiguity. 

The operators 9)o(t) ( OC| ) and hk{t) (4.11) depend on the mean Z(t, H) and moments A a (t, ft), i.e., on the 
solution y( N \t,h) of the Hamilton-Ehrenfest equations (3.4). The solutions of equation ( f4.8| ) in turn depend 
implicitly on y( N \t, K): 

$(x,t,h)=$(x,t,h,y( N \t,K)). 

Below the function arguments y( N ^(t,h) or Wo can be omitted if this does not give rise to ambiguity. For 
example, we may put S)o(t) — £>o(i, ^o)- 

In accordance with the expansion ( 2.1 2| ) and (2.11), the solution of equation (4.8) can be represented in the 
form 



$(x , t, h, * ) = n k ' 2 ^ (x , t, h, * ) + o(h ( - N + 1 ^ 2 ), 



(4.12) 



fe=0 



where 



& k \x,t,h, * )eC t h (Z(t, h),S(t,h)). 



In view of ( pT4| ), for the operators {-ihd t + £o(*,*o)} pTTo| ) and h (k+2 ^ 2 fi k (t, * ), k = 1,N ( pOl] ) the 
following is valid: 



h (k + 2)/2^ Wq) . -pt h (z(t, h),S(t, h), i) -> v f h (z(t, h),s{t, h),i + k + 2), 
{-ihd t + h (t, *o)} : vl(z(t, h),S(t, h), i) -> vl(z(t, h), s(t, h),i + 2). 



(4.13) 



Substitute ( 4.12| ) into (|4.8|) and equate the terms having the same order in ft 1 / 2 , h — > in the sense of ( |OJ ). 
As a result we obtain a set of recurrent associated linear equations of order k to determine the functions 
$«(f,i,ft,*o): 



h 1 {-ihd t + ho(t,^ )W a) =0, 

{- ! M t+ ^(( i $ )}d,W + ^i(i,$o)$ (0) =0, 

h 2 {-ihdt + ho(t, * )}$ (2) + hhi(t, * )* (1) + h 3/2 h2(t, * )$ (0) = o, 



(4.14) 
(4.15) 
(4.16) 



It is natural to call equation ( 4.14 ) for the principal term of the asymptotical solution as the Hartree type 
equation in the trajectory-coherent approximation in mod ft 3 / 2 . This equation is the Schrodinger equation with 
the Hamiltonian quadric with respect to the operators p and x. 
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5 The trajectory-coherent solutions of the Hartree type equation 



The solution of the Schrodinger equation with a quadric Hamiltonian is well known [ £4| p5| . For our purposes, 
it is convenient to take quasi-classical trajectory-coherent states (TCS's) |21j] as a basis of solutions to equation 
(4.14). We shall refer to the solution of the nonlinear Hartree type equation, which coincides with the TCS 
at the time zero, as a trajectory-coherent solution of the Hartree type equation. Now we pass to constructing 
solutions like this. 

Let us write the symmetry operators a(t, ^o) of equation (4.14), linear with respect to the operators Az , in 
the form 

o(t,¥o)=JV«<&(t,*o),AJ}, (5.1) 
where N a is a constant and b(t) is a 2n-space vector. From the equation 

da{t) 



ih- 



dt 



[8o(i,¥ ), «(*)]- =0, 



(5.2) 



which determines the operators a(t), in view of the explicit form of the operator $)o(t, ^o) (4. 1C ) , we obtain 

-ih(b(t),Az) + ih(b(t), Z(t, h)) + 
+ [{-S(t,h) + (P(t,h),X(t,h)) + (Z(t,h),JAz) + i(Az,io 2Z (i,*o)Az)}, (b(t),Az) 

Taking into account the commutative relations 

[Azj, Azk] = ihJjk, j, k = 1, 2n, 



= 0. 



which follow from (1.4) and (1.5), we find 

-ih(b(t), Az) + ih(Az, $) zz (t)Jb(t)) = 0. 

Hence, we have 

b = $) zz (t,y )Jb. 

Denote b(t) = —Ja(t). Then we obtain for the determination of the 2n-space vector a(t) from (5.3) 

a = J9) zz (t,^ )a, ■ 



We call the set of equations (5.4), by analogy with the linear case |18|, a set of equations in variations. 
Thus, the operator 

a(t) = a(t, 9 ) = N a (b(t), Az) = N a (a(t), JAz) 



(5.3) 
(5.4) 

(5.5) 



i s a s ymmetry operator for equation ( 4.14 ) if the vector a(t) = a(t, ^o) is a solution of the equations in variations 

For each given solution Z(t, h) of the Hamilton-Eh renf est equations (3.4), we can find 2n linearly independent 
solutions ak{t) G C 2n to the equations in variations ( |5.4| ). Since to each 2n-space vector au{t) corresponds an 
operator fifc(t, ^o), we obtain In operators n of which commutate with one another and form a complete set of 
symmetry operators for equation (4.14). 

Now we turn to constructing the basis of solutions to equation (4.14) with the help of the operators &k(t, ^o)- 
Equation (4.14) is a (linear) Schrodinger equation with a quadric Hamiltonian and admits solutions in the form 
of Gaussian wave packets 



#o) = N n exp{^S(t,K) + iMt) + iH4>i(t) + (P(t,h), 



Ai 



1 



(Ad 



,Q(t)A£}]}, 



(5.6) 



where the real phase S(t,ft) is defined in ( |4.7| ), Nn is a normalized constant, and the real functions <fio(t) and 
<fii(t) and the complex n x n matrix Q(t) are to be determined. 



Remark 5.1 Asymptotical solutions in the form of Gaussian packets ( |5.6| ) for equations with an integral nonlin- 
earity of more general form than (1.1) were constructed in In this case, the Hamilton-Ehrenfest equations 
depend substantially on the initial condition for the starting nonlinear equation. 
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Substitution of (j5Jj) into ( |4.14| ) yields 

i>{s(t, ft) + z0 o (*) + iWi(t) + (P(t, K), Ax) - (P(t, h), X(t, h)) + -{Ax, Q(t)Ax) - 

-(Ax, Q{t)X(t, h)) - S(t, h) + (P(t, h),x(t, h)) + (X(t, h), Q(t)Ax) - (P(t, h), Ax) + 
+-{(Ax,fj xx (t,y )Ax) + (Ax,S) px {t,q> )Q(t)Ax) + 

+ ([-ihV + Q(t)Ax],S) px (t, V )Ax) + ([^hV + Q(t)A%,Si„(t,* )[-ihV + Q(t)Ax\)}} = 0. 
Equating the coefficients at the terms with the same powers of the parameter h and the operator Ax, we obtain 



(Axfh 1 

(Ax^n 
{Axfh 



iMt) 



Sp [Sj px {t, *„) + ftpp(t, *o)Q(*)] = 0; 



(Ax,0) = 0; 

(Ax, [Q{t) + Si sx (t, * ) + Hxp(t, *o)OW + Q(t)?) px (t, * ) 
+Q(t)Sj pp (t,^o)Q(t)]A£) -0. 



As a result we have 



Sp[Sj px {t) + S}pp(t)Q(t)]dt. 



The matrix Q(t) is determined from the Riccati type equation 

Q(t) + fi xx (t) + Q(m px (t) + fi xp (t)Q(t) + Q(m pp (t)Q(t) = 0. 



(5.7) 
(5.8) 

(5.9) 



Thus, the construction of a solution to equation ( 4.14 ) in the form of the Gaussian packet (5-6) is reduced to 
solving the set of ordinary differential equations ( |5.9|) . 

Let us now construct the Fock basis of solutions to the (linear) Hartree equation in the trajectory-coherent 
approximation (4.14). This is the first step in constructing the solution to recurrent equations (4.14 )— ( 4.16). 

Consider the properties of the symmetry operators Sfe(f) (5.5) of the zero-order associated Schrodinger 
equation (4.14), which are necessary to construct the Fock basis. 

Statement 5.1 Let eti(i) and a 2 (t) be two solutions of the equations in variations and a\(t) and a 2 (t) be the 
respective symmetry operators of equation ( 4.14[ ), defined in ( |5.5| ). Then the equality 

[ai(t),a 2 (t)} = ihN 1 N 2 {a 1 (t),a 2 {t)} = ihN 1 N 2 {ai(0),a 2 (0)} (5.10) 

is valid. 

Actually, upon direct checking we are convinced that 

[a 1 (t),a 2 (t)]=N 1 N 2 [(a 1 (t),JAz),(a 2 (t),JAz)) = 
= ihN x N 2 (Jai(t), JJa 2 (t)) = 
= ihNxN 2 (ai(t), JJJa 2 (t)) = ihN 1 N 2 (a 1 (t), Ja 2 (t)) = 
= ihN 1 N 2 {a 1 (t),a 2 (t)}. 

Here, we have used the rules of commutation for the operators Az. The skew scalar product holds and, hence, 
the statement is proved. 
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Remark 5.2 If the initial conditions for the equations in variations are taken such that 

{ aj (0),o fc (0)} = {a*(0),<(0)} = 0, {a,(0)X(0)} - id k S kj , 
d k > 0, k, j = 



(5.11) 



and Nk = 1/y/hdk, then the following canonical commutation relations for the boson operators of "creation" 
(ctu{t)) and "annihilation" (dfc(i)) are valid: 



Mt),aj{t)] = [at (t), a}(t)], [a k (t),a+(t)} = 6 kj . 



The simplest example of initial data satisfying the conditions (5-11) is 



oi(o) = (6i,o,..., o,i,o,...); 

o 2 (0) = (0,6a,..., 0,0,1,...); 



(5.12) 



(5.13) 



(5.14) 



(5.15) 



Here, we have dk = 21m bk > 0, k = l,n. 
Theorem 5.1 TTie function 

x exp |1 [tf(t, ft) + fi), Af) + i(Af, Q(i)Af)] }, 
where N% = [(Trh)~ n det Do] 1 ^ 4 & a "vacuum" state for the operators cij(t), such that 

aj (t)\o,t) = 0, j = T~rI 

Proof. Actually, substituting ( |5.5| ) and ( 5.14 ) into ( |5.15| ), we get 

\0,t)[(Z 3 (t),Q(t)Ax) - (#,•(*), Ax)] =0, 

since we have 

Q(t)^(t) = BCtJCr 1 (*)£,-(*) = Wj(t). 

Recollect that from the fact that the matrix Do is positive definite and diagonal follows det C(t) ^ 0, and 
so the matrix Im Q{t) is positive definite as well (see Appendix A). 

Let us define the denumerable set of states \v,t) as a result of the action of the "creation" operators upon 
the "vacuum" state |0, t): 

n - 

|M) - M,*o) - -r(a+(i,tf )r|0,t,*o) = TT _(a+(t,* o )) I/f '|0,t,*o). (5.16) 
iA - LJ - i/j.! 



fc=l 



By analogy with the linear theory (k — 0), we call the functions \u, t) ( |5.16 ) quasi- classical trajectory-coherent 
states and consider their simplest properties. 



Statement 5.2 The relationships 



a,k\v,t) = y/vk\v[ \t), 
a+\v,t) =V^k~TT\v ( k + \t), 

= (^1 ± ^1,*> ^2 ± $2,k, ...,V n ± 5 n<k ) 



(5.17) 



are valid. 
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Actually, we have 
It follows that 



n 1 n 

a>,i) = ]J [a j7 (a+r}\0,t) = ]J^(atr- S ^\0,t) = 



k=l 



'V k \ 



k=l 



TT 



(a+)^-^|o,t) = ^l^ ) ,t); 



n 1 



and thus the statement is proved. 

Statement 5.3 The states \v,t, ^o) t£R and S T 5 ^ form a set of orthonormal functions: 

(*o,*Vk*>*o) = <5^<, e Z'j. 

Let us consider the expression 

1 



(5.18) 



(*o,*yif,*,*o> = 



vVR 



(*o,i,0|a u (t,* )[5 + ]' y (i,*o)|0,t,*o)- 



Commuting the operators of "creation" and "annihilation" in view of commutation relations (5.12) and using 
relationship (5.15), we obtain 

{t,i/\v,t) = {t,Q\Q,t)8 u y. (5.19) 

Then we calculate 



(*,0|0,t) 



| det C(t)| 



f 2 

/ exp[— — Im S(x, t)]dx. 



(5.20) 



In view of (A. 8) and the explicit form of the complex phase in (5.14), we have 

Im S(x,t) = -(Ax,Im Q(t)Ax), 



The matrix Im Q(t) is real and positive definite; hence, the matrix y'lm Q(t) does exist, such that 



Let us perform in the integral of ( |5.20 ) the change 



|*= 4= Vim Q(t)Ax 



and then obtain 



(*,0|0,*) = 



Vdet D J 



Vdet D 



since detZ?o = II I m °k- Thus, the functions \v,t, ^o) ( 5.16| ) form the Fock basis of solutions to equation 



fc=i 



Q4.14Q , Q.E.D. 

Theorem 5.2 Let £/ie symbols of the operators T~l(i) and V(t, VP) satisfy the conditions of Supposition The 



for any v £ 17\ the function 



® v (x,t,h) = \v,t), 



(5.21) 



where the functions \u,t) are defin ed by formula ( 5.16| ), is an asymptotical (accurate to 0(h 3 / 2 ), h — > 0) solution 
to the Hartree type equation (1.1) with the initial conditions 



v v (x,t,h)\ t=0 = \v,t)\ t = . 



(5.22) 



1G 
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The principal term of the quasi-classical asymptotic of the Hartree 
type equation 



The solution of the Cauchy problem (1.1), ( |5.22 ) is a special case of the quasi-classically concentrated solutions 
of equation (1.1). However, in the case of arbitrary initial conditions (2.18) belonging to the class V\, the 
functions \ v 1 1) are not asymptotical solutions of the Hartree type equation (1.1). This is a fundamental difference 
between the complex germ method for the Hartree type equation ([O]), being developed here, and a similar 
method developed for linear equations |l^, ?, |2l] |. The coefficients of the Hartree type equation in the trajectory- 
coherent approximation (4.14) depend on the initial condition fl2.18| ) since they are determined by the solutions 
of the set of Hamilton-Ehrenfest equations. It follows that among the whole set of solutions to equation (4.14) 
only one (satisfying the condition ty(x,t, h)\ t= o = ^o(3,h)) will be an asymptotical (accurate to 0{h A/z )) 
solution to equation (1.1). However, the Fock basis (TCS's) \i/, t) makes it possible to construct an asymptotical 
solution to equation (1.1) with any degree of accuracy in ft 1 / 2 , h — ► in an explicit form which would satisfy 
the initial condition ( 2.18 ). 

Let us illustrate in more detail the relation of the solutions of the associated linear Schrodinger equation to 
the solution of the Hartree type equation. To do this, we construct the Green function of the Cauchy problem 
for the zero-order associated Schrodinger equation. Although the Green function G^ (x, y, t, s) for quadric 
quantum systems is well known [^5[ [38] [39], |l|, we give for completeness its explicit form, as convenient to us. 
This function will allow us to demonstrate explicitly the nontrivial dependence of the evolution operator of the 
associated linear equation on the initial conditions for the starting Hartree type equation. 

By definition we have 

[-ihdt + ho(t *o)]G(°) {x, y, t, s, tf o) = 0, 

limG(°)(x,£M,#o) = 8(2 - y), (( '"" 



where the operator F>q is defined in ( 4.10 ). We make use of the simplifying assumption that 



det Sjpp(s) ^ 0, det 



dp k (t,z ) 



7^0. 



(6.2) 



If the condition (f3j) is not valid, the solution of the problem can be found follow ing the work j39] [to) . For 
the problem under consideration, exact solutions of the Schrodinger equation (4.14) are known: these are the 
functions \v, t, \&o) ( 5. 16| ) that form a complete set of functions. Thus we have 



E 

IH=o 



(6.3) 



where 



® v (x,t,h) = M,*o)- 



Details of similar calculations can be found, for instance, in p5[ . However, for our purposes the following 
approach seems to be convenient. 

Let us carry out an ft -1 Fourier transform in equation d6.1|). For the Fourier transform of the Green function 



GW(p,y,t,s,V ) 



d.f 



^G (0) (x , y , t, s, * ) exp { - % - (x, pj } 



we obtain 



ihd t + hoik 3, t, *o)]G(°> (p, y, t, s, o) = 0, 
1 



limG(p,y,t,s^ ) = X27Tm)n/2 



exp 



(6.4) 



(6.5) 



Here, p — p and x — ihjL and the symbols of the operators of equations (3.1) and (6.5) coincide: 



fto{p,x,t) = Sjo(p,x,t). 
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Equation (5.5) coincides to notations with ( 4.14 ) and, hence, admits solutions of type (5.6) 

(p, y, t, s, * ) = exp { - % - [s (t, s, y) + {G(t, s,g),Ap) + - (Ap, Q(t, s, y)Ap)} } , 

where Ap — p— P(t, h). Here, the functions So(t, s, y) = So(t), G(t, s, y) = G{t) and Q(t, s, y) = Q(t) 
determined and, according to (|6.5|), satisfy the initial conditions 



lim Q(t, s, y) = 0, lim G(t, s, y) = y, lim S (t, s, y) = (p , y)- 

t—>s t — >s t — >s 



Substituting ( |6.6[ ) into (6.5), we write 

&°\p,y,t,s,* ){ - S (t) - (d(t),Ap(t)) + (G(t),P(t,h)) - ±(Ap,Q(t)Ap)+ 
+{P(t, ft), Q(t)Ap) - S(t, h) - (P{t, ft), (G(t) + Q{t)Ap- P(t, ft))) + (X(t, ft), Aj3) 
+\ [((G(t) + Q(t)Ap~ P(t, h)),f) xx (t)(G(t) + Q(t)Ap- P(t, ft))) + 
+((G(t) + Q(t)Ap- P(t,K)),f, xp (t)Ap) + (Ap,f) pp (t)Ap) + 
+(Ap,9) px (t)(G(t) + Q(t)Ap- P(t,h)))] +yS P [Sj xx (t)Q(t) + fi xp (t)]} = 0. 

Equating the terms with the same powers of Ap, we obtain the following set of equations: 

-Q = ¥) px {t)Q + QSj xp (t) + Qft xx {t)Q + S) pp {t) = 0, 

-{t - P{t, ft)) + Q(t)S) xx {t)(G ~ P{t, ft)) + Sj px (t)(G - P{t, ft)) = 0, 

-S + (X(t, h),P(t, ft)) - S(t, ft) + y Sp [$j xx (t)Q(t) + 9) xp (t)]+ 

+\0(t) - P(t, h)),f) xx (t)(G(t) - P(t, ft))) = o 



with the initial conditions (6.7) 



Let B(t) and C(t) be solutions of equations in variations (A. 2) with the initial conditions 

B(t)\ t = s = B (s), C{t)\ t=s = 0, B t (s) = B (s) 

and the matrix Im Bq(s) be positive definite. 

In view of ( |6.2| ), the solution to the Cauchy problem (A. 2), ( |6.11 ) will then have the form 



B{t) = X\(At)B (s), C{t) = -\\{At)B {s), At = t-s, 



where the matrices A| (t) and \\(t) are defined in (A. 16). The matrix 



Q(t) = C{t)B-\t) = -\\{At){\l\At)Y 



will then satisfy equation ( J6.8| ) with the initial conditions (6.7). 
Provided that (|J) is valid, from ( |A.12j ) and ( |A.10D follows 



G{t) = {B-\t)Y (y - x ) + X(t, ft) = X^(At)(y- x ) + X(t, ft). 
In a similar manner, we obtain for So 

t 

S (t,s,h) = S(t,h)- S(s,ft) + y j drSp [S) xp (t) + Sj xx (t)Q(t)}+ 



(6.6) 
are to be 

(6.7) 



(6.8) 
(6.9) 

(6.10) 

(6.11) 

(6.12) 
(6.13) 

(6.14) 



+ 2 / dr((G(r) - X(T,h))^ xx (T)(G(T) ~ X(r,h))) + (p Q ,y). 



(6.15) 
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In view of (A. 12) and Liouville's lemma A. 3, we obtain 



dr Sp [Sj xp {r)+Sj xx (T)Q(T)} = -lndetB" 1 ^)!* = 



= l/detB (s)^ 



/ det tf {s) \ 
V detB(t) l 



In det A 4 (At). 



To calculate the last integral in (6.16), we use relationship ( A.15| ) and, in view of ( |6.15| ), we get 



dr((G(r) - X(r, ft)), 9j xx (r)(G(r) - X(r, ft))) 
= ^{(y-S ),X 2 (At)X^ 1 (At)(y-x )}, 



where the matrix A| (t) is defined in (A. 16). Hence, we have 



(6.16) 



(6.17) 



S (t,s,h) = S(t,h)-S(s,h)-—liL(det\i(At))+ 

+ ^({y - x ), \2{At)\- 1 (At)(y - x )) + (po,y)- 
Substituting §1$) , (|1J), and (|o|) into (|J), we obtain the well-known expression (see, e.g., |39[| ) 



: exp 



(- l(S(t,h) -s(s,h)- 
i ft 



{2wih) n / 2 v /detA 4 (At) 
2 (At) A J 1 (At) (y — x )) 

-^(Ap,X^(At)(y-x )) + ^(Ap,A 4 - 1 (At)A 3 (At)Ap)}. 



^{{y-x ),X 2 (At)X^(At)(y-x )) - % -(j> Q ,y) - '-(Ap, X(t, ft))- 



Now we substitute (6.19) into (6.4) and make use of the relationship 



dx exp 



(x, Fx) + (6, 2>] = V(2^) n detr-! exp { ^ } 



(6.18) 



(6.19) 



(6.20) 



in which we put T = -(i/h)X^ 1 (At)X 3 (At), b = -(i/H)[X(t,h) - x + A~ 1 (At)(y ~ x )]. We then obtain 

(f , y, t, s,* )= 7= = = feX , == exp { I [,S(t, ft) - 5(.s, ft)+ 
A/det(-^27rftA3(At)) LftL 

+(P(t, ft), Ax) - (p , (y- Sb)> - 3d), A x (At)A3 x (At)(? - f )>- 



--(A^Aa^At)^-^)) - 2< A ^ A,^ 1 (At)A 4 (At)Af)]}. 
Here, we used the relationships 

A<(t)A 4 (t) - A 3 (t)A 2 (t) = I„ xn , A 3 (t)A*(t) - Mt)K(t) = 0, 



(6.21) 



that follow immediately from ( A. 5 ), ( A. 19 ) , and from the definition of matriciant ( A. 16 ). 
Let us consider the limit of expression ( |6.21[ ) for At = t - s ->• 0. We obtain 

Ai (At) = l nxn + O(At), A* (At) = -Sj pp (s)At + 0{{At) 2 ), 



(*) + 0((At)°), A 4 (At) = I„ x „ + O(At), 
A 2 (At) = O(At). 
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It follows that for short times we have (see, e.g., jyj) 

lim G (0) (s, V,t, s,^! ) = , 1 =x 

x exp{^ Ki ((x-y),Sj-i(s)(x-y)) + 0(At°)}. (6.22) 

Thus we have proved 

Theorem 6.1 Let the symbols of the operators Ti.(t) and V(t, "P) satisfy the conditions of Supposition ^J. Then 
the function 

¥°\x 1 t 1 h) = U^(t,0,^ o )^ o , (6.23) 



where (t, 0, ^o) is the evolution operator of the zero-order associated Schrodinger equation ( 4.14 ) with the 
kern el G^(x,y,t,0, Vl'o), is an asymptotical (accurate to 0(h 3 ^ 2 ), h — > 0) solution of the Hartree type equation 
(1.1) and satisfies the initial condition 

*< >(j?,t,ft)| t=0 = *o. 



Remark 6.1 TTie principal term of the quasi-classical asymptotic ^>^°\x,t,h) will not change (accurate to 
0(h 3 / 2 ), h — ► 0) «/ i/ie phase function S^ N \t, K) in the operator S)o(t, ^o) is substituted by its value S 2 (t, H) for 
N = 2 and we restrict ourselves to the first terms in h — > in i/ie phase trajectory Z^ 2 '(t,h), and in the other 
expressions Z( N \t,h) is changed by Z°(t,h). 



7 Quasi-classically concentrated solutions of the Hartree type equa- 
tion 



Now we construct asymptotical solutions to the Hartree type equation (1.1) with an arbitrary accurac y in 
powers of Vh. To do this, we find asymptotical solutions to the associated linear Schrodinger equation ( |4.6| ) 
with an arbitrary accuracy in powers of \fh. Let us present an arbitrary initial condition $o(2?i K) € V2 



N 



<p (x^) = z 2 hk/2 4 k) &fr), 



(7.1) 



fe=0 



where 



ft) eC£(z„, So). 

Then for the recurrent associated linear equations (4.14 )— ( 4. 16| ) we arrive at a Cauchy problem with initial data: 

$ (fc) | t=0 = <s> { o ] (x,h) fc = rvv. 

The solution to these recurrent equations can readily be constructed as its expansion over the complete set of 
orthonormalized Fock functions \u,t) (5.15). As a result we obtain 

oo 

^°\x,t,h) = M,*o)<*o,0,H$ o O) (f,ft)), (7.2) 



\u\=a 

OO - 

% - |M,*o) / dr(vl/ ,T,H^i(t,*o)<i> (0) (f,r,ft)), 
l-l=o o 

OO 

$W(x,t,h) = k,t,*o)(*o,o,H* 2) (^fi))- 

\v\=0 



(7.3) 
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oc „ 

-\ E l^*'*o) / dT<* ,r,i/|fli(t,*o)^ (1) (^r,ft))- 
|H=o o 

OO - 

~\ Yl M>*o) / dr(^ ,r,^2(t,*o)$ (0) (f,r,?i)), (7.4) 

Denote by ^ N ^(t, ^o) the operator denned by the relationship 

t 

£W(t,¥ )*(*,i)= /<fr#o(*,T,* )tf W fo r), (7.5) 



n 



where J7o(i, t) is the evolution operator of the associated Schrodinger equation ( 4.14 ) and the following notation 
has been made: 

AT 

h {N) (t,y Q ) = ^h k / 2 hk{t,y ). (7.6) 

Thus we have proved the following statement: 

Theorem 7.1 Let the symbols of the operators Ti.(t) and V(t, \P) satisfy the conditions of Supposition T/ien 
i/ie function 

n . k 

*W(*,t,ft) = £_{:j^W(t,$ )} £> (t,0,*o)*o(^fi)., (7.7) 

fc=0 

where N > 2, is an asymptotical, accurate to 0(h^ N+1 ^ 2 ), solution of equation (|l.l| ) and satisfies the initial 
condition ( |2.§| ). 

8 The Green function and the nonlinear superposition principle 

Let us show that in the class of trajectory-concentrated functions for the Hartree type equation ([hi]) we can 
construct, with any given accuracy in h 1 / 2 , the kernel of the evolution operator or the Green function of the 



Cauchy problem for equation (LI). The explicit form of the quasi-classical asymptotics ^ N '(x,t,h) (7.7) 



makes it possible to obtain an expression for the Green function G^- N \x, y, t, s, ^o) valid on finite time intervals 



t £ [0, T]. Actually, according to (7.7), for any function ip(x,h) € the solution of the Cauchy problem with 
the initial condition 

3>(x,t,h)\ t=0 = ip(x,K) (8.1) 



for the associated linear Schrodinger equation fl4.8|) has the form 

& N \x,t,h) = 



flW(t,* )/ dyG^(x,y,t,0,^My,h) + Oih^/ 2 ) = 
djJGW (x, y, t, 0, *„Mi7, h) + 0(^ A '+ 1 )/ 2 ), 



where 



N 



&*{t,v ) = ^±{±j<*(t,*o)} ' (8 - 2) 

fc=0 

and the function G^ (x, y, t, s, ^>o) is defined in ( |6.2l| ). 
It follows that 

(x, £ i, 0, ¥«,) = R {N) (t, *o)G (0) (2, y, t, 0, *„) . 
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Since R^ N \0, 'Jo) = 1, we have for an arbitrary s < t 

G^(x, y, t, s, , * ) = R (N) (t, * )G(°»(f, y, t, s, *o)(R (N) (s, *o)) + 



(8.3) 



being the Green function of the Cauchy problem (8.1) with s ^ 0. Obviously, for the functions G^ N \x, y, t, s, ^o) 
the following composition rule is valid: 



duG^\x 1 u 1 t,T,^o)G^\u,y,T, s ^ ) = G ( - N \x,y,t,s,^ )+ 

+ 0(^ w+1 '/ 2 ). 



(8.4) 



Denoting by U^ N '(t,0, V^o) the approximate evolution operator of the linear equation (4.; 

uW(t,o,v Mx,h)= [ dyGW(x,y,t,0,VoMy,h), 



we obtain it from (8.3) in the form of the T-ordered Dyson expansion 



N ■ k r 

^W(i,0,* ) = $3(-^) / d fc ^i(T L ,t J *o)-"«i(n,*,*o)t'b(<,0 > *o). 



fe=0 



(8.5) 



Here, we have used the following notations the domain of integration is an open hypertriangle 

A> = {r S [0, t] k \t > n> t 2 > ■■■ > t n > s}; 

the operator f)i{r,t, ^o) is a perturbation operator in the representation of the interaction 

»i(r, t) - #<,(*, r, * )^ W (r, *o)^ + (r, t, * ), 



(8.6) 



where f)( N ' (t^o) has been defined in ( |4.9|) , and E/pft, s , ^q) is the evolution operator of the associated linear 
Schrodinger equation wit h th e kernel G^Hx, y,t,s, ^ p) (|6.2l|) . 

In view of Statement 4.1, the action of operator ( |3.5] ) on the function (p = ^q (x, ft) determines the asymp- 
totical solution of the Cauchy problem (fTl|)-( |2~8| ) for the Hartree type equation (1.1) 



*W(x,t,h) = #W(t,0,tf )*o(2,fi), *o(2,») G V° h . 



•7) 



It follows that operator ( p.5[ ) is an approximate evolution operator for the Hartree type equation (1.1) in the 
class of trajectory-concentrated functions. 

For the constructed asymptotical solutions, from expression ( |B.7| ) immediately follows |45], [l6j 

Theorem 8.1 (nonlinear superposition principle) Let \& i (of, t, ft, y[ N \t, ft)) be an asymptotical, accurate 
to 0(H( N+1 ^ 2 ), solution of the Cauchy problem for the Hartree type equation ( |l.l| ) with the initial condition 
^oi(x,h) an d the function ^(^j t, ft, (t, ft)) is a solution of the same problem with the initial condition 
^02(x,h). Then the solution of the Cauchy problem with the initial condition 



y 03 (x,h) = C 1 ^ 01 (x,h) + C 2 ^Q2(x,h), 



c\,C2 — const, 



has the form 



y 3 (x,t,h,yi N \t,h))=m N ~>(t,0,y 03 )y 03 (x) = 

= Ci m N \t, 0, *03)*0l(£) + C2U {N) (t, 0, *03)*02(aO = 

= c x *i(f , t, h, y { 3 N \t, ft)) + c 2 ^ 2 (x, t, ft, yi N) (t, ft)). 

Here, 2/1 (i,ft) denotes the solution of the Hamilton-Ehrenfest equations of order N, N > 2 (3.4) with an 
initial condition which is determined from the functions ^ok(x,h), k — 1,3, respectively. 
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9 The one-dimensional Hartree type equation with a Gaussian po- 
tential 

Let us illustrate the above scheme for constructing asymptotical solutions by the example of a nonlinear in- 
teraction with a Gaussian potential |l7|. By this example we shall show in an explicit form how the proce- 
dure of constructing quasi-classically concentrated solutions to the Hartree type equation necessarily leads to 
Hamilton-Ehrenfest equations. Moreover, it becomes possible to elucidate the "nonlincarity" of the generalized 
superposition principle for the Hartree type equation. 

We write equation (1.1) with a Gaussian potential for the one-dimensional case as 

■{x-yf- 



VP) 

-ihd t + — h >cVq I dy exp 

2m 



2 7 2 



\V(y,t)\ 2 }V = 0. 



(9.1) 



In this case, for the class of functions Vt(S(t,H),Z(t)) in which we shall find solutions to equation (9.1), in 
accordance with (O) , we find 



Vl{S(t,h),Z(t)) = 

* : $>(x,t,h) =<p(-j!L,t,h)exp Us(t)+P(t)Ax) 



(9.2) 



Ax 



Here, the function ip(£, t,K) G § (Schwartz space) with respect to the variable £ = — — and depends regularly 

on h with Ax = x — x(t). The functions S(t,H) and Z(t) = (P(t), X(t)) are real, depend regularly on H, and 
are to be determined. 

Let us expand the exponential in equation ( |9.1[ ) in a Taylor series for Ax = x — X(t), Ay = y — X(t) and 
restrict o ursel ves to the terms of the order of not above four in Ax and Ay. In view of the estimates (p-7|), 
equation (9.1) will then take the form 



-ihdt + 



P 2 {t) P{t)Ap Ap 2 



2m 



2m 



-xVn 



1 - ^(Ax 2 - 2Ax«i 1} (t, K) + af (t, ft)) 



2 7 2 



--^(Ax 4 - 4Ax 3 a4 1} (t, K) + 6Ax 2 a y "> (t, h) 
87 

-AAxaf{t,h) +a ( £ ) (t,h))] W = 0(H 5 / 2 ), 



(9.3) 



where Ap = p — P(t), h = x||* || 2 , and 

a * ) (*.^)=|i J {Ay) k \^{y,t)\ 2 dy, k = 0^ 



+00 



are the fc-order moments centered about X{t). Equation (3.3) can be simplified if we make the change 



^(x, t, K) — exp 



P 2 {t) 



2m 



%V + -^V a xx [t,h) dt\<f>(x,t,h), 



(9.4) 



where 



+00 



(2) 

o xx (t,K) = a\,'{t,h) = 



1*1 



Ay 2 \*(y,t)\ 2 dy 
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is the variance. The function <j>(x, t, K) belongs to the class V\\S(t, H),Z(t)) , where 



S(t, h) = S(t, K) 



P 2 (t) 



2m 



+ xV - -^Voa xx (t, h) 



and satisfies the equation 



P(t)Ap Af 
-iWt H h h xVb 

TO ZTO 



^(A a ; 2 -2Axai 1) (t,ft)) 



--i-(Ax 4 - 4Ax 3 a4 1) (t, h) + 6Ax 2 a4 2) (t, R) 
87 

-AAxa ( i\t, h) + h)) ^$ = 0{tf' 2 ). 



Here, we have made use of 



(9.5) 



a$\t,h)=o£ ) (jt,h), k = l,N. 
We shall seek the approximate (mod/1 5 / 2 ) solution <J>(x, t, ft) to equation ( |9.5[ ) in the form 

t, ft) = $ (0) (x, t) + \/ft$ (1) {x, t) + ft$ (2) (x,t) + ..., 

where $( fc )(ir,i) g C^(S(t,h),Z(t)). In equation ( |9.5| ) we equate the terms having the same estimate in \/ft 111 
the sense of (2.7). Denote by Lq the operator 



L = -ihd t + -P(t)Ap + ^-Ap 2 - ^| A* 2 ; 

to 2to 27^ 



Earlier we have shown that Lq = 0(h). We then have 



ft 1 : (L + ^Axa%)^=0; 



ft 3 / 2 : v^io + ^Axa^))^ 1 ) = -2Vft^AxRc($W|Ax|<i>«)$ 



(0) 



2 7 2 



ft(L + ^A^)*! 2 ) = -2ft^Ax 
+ ($ (1) |Aa;|$ (1) )l$ (0) +2Re($ (0) |Ax|$ (1) )$ (1) U$(t)4 
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2Re($^|Ax|$( 2 ) 



s(0) 



The function 



1/2 



i /irt 



1 



exp ^l^ + wAx+ 2 Aa;2w c(t)y 



\\C(t)\ 



exp 



ftl m 2 CW 



(9.6) 



is a solution of equation ( |9-6| ). Here, we have used the fact that X(t) and P(i) are solutions of the ordinary 
differential equations 

f p = 



I x=— , p(0) = p ,#(Q) = so, 

k TO 

and C(i) denotes the complex function satisfying the equations 



(9.7) 



C = 



B(t) 



(9.8) 
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Equations 



are Hamiltonian equations with the Hamiltonian TC(p,x,t) — p 2 /(2m) and their solution is 



P(t) = mv = po, X(t) = — t + xq 



Similarly, equations (9.7) arc Hamiltonian equations for a harmonic oscillator with frequency 



' &\Vo\ 
rwy 2 



and their solution is 



C(t) — c\ exp ( 
B(t) = mC(t). 

For the initial conditions ( 5.1 3| ) 



x\Vo\ 
rwy 2 



t) + c 2 exp ( 



*\Vo\ 

mrf 1 



two cases are possible: 



C(0) = 1, B{0) = 6, Im b > 



1) kV o >0: C(t) = ch(fit) + -sh(fii), 



2) xV < : C(t) = cos(Ot) + - sin(fii). 



(9.9) 



In the linear case (x = 0), the frequency 17 = and equations ( J9 .8| ) become equations in variations for equations 
(9.7). In view of (|9.8|), we find the variance of the coordinate x in an explicit form: 



cr xx (t, K) 



; Im b 



71 H 



\C(t)\ 



exp 



m , n Im & 
At 

ft |C(*)I 2 



da; = 



\C(t)\ 2 h ^ 
2ro • Im 6 



(9.10) 



Then we get 



exp 



2 



2m ■ Im 6 



|C(i)| 2 cft) *<°> (*,*>&) 



It can readily be noticed that a)^' 0) (t, K) = a^ (0) (t, K). Hence, from ( p.S| ) and (|9. 10[ ) it can be inferred that for 

(2) 

kVq < the variance o>\^ 0) (t,H) is limited in t, i.e., we have |a" ax (t)| < M, M = const, and for i<Vb > it 
increases exponentially. In the limit of 7 — > and with Vq = (2^7) _1 / 2 , equation (9J3) becomes a nonlinear 
Schrodinger equation, while in the case where hVq < (>cVq > 0) it corresponds to the condition of existence 
(nonexistence) of solitons. Note that if a^( 0) (t, H) = 0, the equation for the function takes the form 



Ln ■ $(°) = 



(9.11) 



becoming a Schrodinger equation with a quadric Hamiltonian. We shall find the solution to equation (9.11) 
satisfying an additional condition a^ 0) (t,h) = 0. To do this, we denote 

a(t) = N a (C(t)Ap - B(t)Ax). 



If C(t) and B(t) are solutions of equations (p.q), the operator a(t) commutates with the operator Lq. So the 
function 

*L 0) = li^m^ 
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will also be a solution of the Schrodinger equation (9.5). Commuting the operators a + (t) with the function 
$q (x, t, ft), we obtain the Fock basis of solutions for linear equation ( p. 11 ) 



^(x,t) = ±N k ^(x,t)(-i) k [C*(t)] k 



. dx 



2m Im b 

MP 



-Ax 



1 



±N^\x,t)(-i) k [C*(t)] k 



( \J hm Im b \ k ( . V m Im b \ 



V |C(*)| 



|c(t)i^; ! 



where H n (£) are Hermite polinomials. Determining N a from the condition [a(i),a + (i)] = 1 and representing 
the solution of the equations in variations as 



C(i) = |C(t)|exp{-iArg [C(t)]}, 



we get 



< ) (^0 = ^H) fc exp{- l fcArg[C(t)]}(i=) i^Ax^^) ^(x s t). (9.12) 



Using the properties of Hermite polinomials, we can readily be convinced that the mean a (0) (t, ft) = 0, k = 0, oo. 
Then we have 



exp 



I ft 



.2 

V2m 



r ill 4- xF )t + xFo^flAar 2 !^)] }$i 0) (x,t, ft). 



Similarly, we find 



(^lA, 2 !^))^^-^ / A, 2 |<I^(M)| 2 ^ 



ft|Cft)| 2 (2fc + l) 
2m Im 



and for functions (t, ft) obtain 



^f\x,t,h) = exp 



+kV 



h(2k + i) 

2m Im b 



\2m / 
|C(r)| 2 dr)]$(°)(x,t,ft). 



(9.13) 



The functions ( 9.1 3| ) are approximate, accurate to 0(ft 3 / 2 ), solutions of the Hartree type equation (9.1). How- 
ever, since for the linear combination 



$(x,t) =C 1 ^ ) {x,t)+C2^\ V, (x,t) 



(0), 



(9.14) 



the condition ct£\t, ft) = is not fulfilled, &(x,t) is not a solution of equation (3.6) and, hence, the linear super- 
position principle is invalid for the functions ( p. 13 ) even in the class of asymptotical solutions V^(S(t, ft), P(t), X(t)) 

accurate t o Q(ft 3 / 2 ). Thus, the presence of the term a^ 0) {t, ft) in equation ( |9.6| ) violates the linear superposition 
principle ( 9.14 ). 

We seek the solution to equation (9T ) in the class VySSit, ft), Z(t, ft)) , i.e., localize the solution asymptotically 
i ii th e neighborhood of the trajectory z = Z(t, ft) depending explicitly on parameter ft. With that, the estimates 
(2.7) remain valid. Let us take the dependence of Z(t, ft) on the parameter ft — > such that the equation for 
the function (x, t, ft) be linear. For doing this, we subject the functions X(t, ft) and P(t, ft) to the equations 



2xVb (i) . . 1 (3) , 

P = — —^(0) (*: h ) + 72 %(o) (*, ft )> 

x = p/m, 



(9.15) 
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and the functions C(t) and B(t) to the equations 



■ - ' 4 7 2"<E.(0)> 



C = B/to. 

The function i, H) will then satisfy the equation 

L * (0) = 0. 



(9.16) 



(9.17) 



Unlike equations ( |9.6| )-( |9~8| ), equations ( 9.15 )-( j9.1?l ) are dependent. Note that, with in the accu racy under 
consideration, the principal term of the asymptotic will not change if equations (3-15) and (9.16) are solved 
accurate to 0(fi 3 / 2 ) and 0{K), respectively. Then equations ( 9.1 5| ) become 



2kV (1 



7 
P 



2 4> 



(9.18) 



and equations (|9.16| ) coincide with equations (|9.8| ) and their solution has the form (p.9j). Equation (|9.17| ) is linear 
and its general solution can be represented as an expansion over a complete set of orthonormalizcd functions 



(9.19) 



k=0 



Here, ^^\x,t) is determined by expression (9.12) where X(t) and P(t) ought to be replaced by X(t,H) and 
P(t,H), respectively. Substitute (p.!9|) into (9.S). In view of the properties of Hermite polinomials 



+ OC 



-\-oo 



_ t 2 ± 



we obtain 



* to lm 6 



X j 



n=0 



VH|C(t)| 

m Im 



.X. ^ 



n=0 



Equations ( 9.18 ) will then take the form 



p = Q l Vh-\C(t)\ 
p 



where 



9l = ^In^^(2 C «+ 1+nC "-0 C ' 1 



TO7 2 Im b 

1 n=0 

Integration of the obtained equations yields 

P{t,h) = mX(t,H), 

t T 

X(t,h)=Vh— [ dr [ \C(s)\ds + ^-t- 

TO J J TO 



Xq. 



(9.20) 



o o 
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As a result, the principal term of the asymptotic can be represented in the form 

P 2 (r,h) 



ex P {-^[/ (p(T,h)X{ T ,h) 



2m 



■ ne- 



|C(r)| 2 )dr]}$(°)(x,t,ft), 



(9.21) 



where 



02 



2m Im b 



oo 



„- -2 + ( n + gK + C™ - U ) C n-2 



It follows that function ( 9.21 ) depends on 0i and O2 as on parameters: 

¥°\x,t,h) = *(°)(x,t,ft,e 1 ,e 2 ), 

Here, 8i and 02 are determined by the sets of equati ons ( |9.20 ) and ( |9.22 ), respectively. 
Let us consider the Cauchy problem for equation (9.1): 



t=o 



t=0 



t=o 



*2o(a:), 



* 3 o(z) = Gi*io(^) + G 2 *2o(a;) 



(9.22) 



(9.23) 



where Gfc = const. Denote by ^(rr, t, ft, 2 ) the principal term of the asymptotic solution to equation 
(9.1), satisfying the initial conditions ( |9.2S ). Then from the explicit form of function ( |9 .13 ) follows 



# 3 (z, t, ft, e?, ei) = Gi* 2 (as, t, ft, el ei) + g 2 * 2 (x, t, ft, e?, e|). 



(9.24) 



Relationship ( 9.24[ ) represents the nonlinear superposition principle for the asymptotical solutions of equation 
(jO) in the class V^(S(t,h),Z(t,h)). 



The work was supported in part by the Russian Foundation for Basic Research (Grant No. 00-01-00087). 



Appendix A. The set of equations in variations 



We already mentioned that to construct solutions to equation (4.14) in the class Vl, it is necessary to find 



solutions to the equations in variations (5.4) and to the type Riccati matrix equation (5.9). Let us show that 
the solutions of the Riccati type matrix equation can be completely expressed through the solutions of the 
equations in variations a(t). 

Let us present the 2n-space vector a(t) in the form 

o(t,* ) = (W{t,* ),Z(t,Vo)), 

where the n-space vector Wit) = Wit, V&o) is the "momentum" part and Z(t) = Z(t, ^q) is the "coordinate" 
part of the solution of the equations in variations. Thus we can write the latter as 



W = -S) xp it, V )W - f) xx (t, * )Z, 



(A.1) 
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The set of equations (A.l) is called a set of equations in variations in vector form. Denote by B(t) and C(t) 
the n x n matrices composed of the "momentum" and "coordinate" parts of the solution of the equations in 
variations: 

B(t) = (# x (t), W 2 (t), . . . , W n (t)), C(t) = (Zt{t), Z 2 (t), . . . , Z n {t)). 
The matrices B(t) and C(t) satisfy the set of equations 



C = S) pp (t,y )B + f) px {t)C, 



(A.2) 



which is called a set of equations in variations (5.4) in matrix form. 

Let us consider some properties of the solutions of this set of equations, which determine the explicit form 
of the asymptotical solution of the Hartree type equation and its approximate evolution operator. 



Remark A.l The set of equations in variations (5.4) is a set of linear Hamiltonian equations with the Hamil- 
tonian function 

H(a,t) = -{a,Sy zz (t)a), ae 



The complex number {at, 02} — {at, Ja 2 ) is called a skew-scalar product of the vectors at and a 2l a^ S C 2 ™. 
Obviously, the skew-scalar product is antisymmetric: 

{at, 02} = -{02, ail- 



Statement A.l The skew-scalar product {at(t) , a 2 (t)} of the two solutions, a\{t) and a2(i), of equations in 
variations (5.4) is invariable in time, i.e., 

{at{t), a 2 {t)} = {ai(0), a 2 (0)} = const , (A.3) 
{at{t), al{t)} = {ai(0), a* 2 (0)} = const . (A.4) 

This statement can be checked immediately by differentiating the skew-scalar product {ai(i), a 2 (t)} with respect 
to t: 

j t {at{t),a 2 (t)} = (di(t), Ja 2 {t)) + (ai(i), Ja 2 {t)) = 

= {JSj zz {t)at{t),Ja 2 (t)) + {at{t),JJSj zz {t)a 2 (t)) = 
= {at(t),fi zz (t)a 2 (t)) - {a 2 {t),Sj zz {t)a 2 {t)) = 0. 



Here, we have made use of the fact that J 2 = — l 2nx2n and J* = — J. Relationship (A4) follows from (A.4) 
since a 2 (t) is also a solution of the equations in variations in view of the fact that these equations are real and 
linear. 

For the set of equations in variations in matrix form, Statement [A.l| becomes 



Statement A.2 The matrices 



D = ±[C + (t)B(t)-B+(t)C(t)}, (A.5) 
D Q = C t {t)B{t)- B*(t)C(t), (A.6) 



where the matrices B(t) and C(t) are arbitrary solutions of the set of equations in variations (A.2), are invariable 
in time, and so we have 

D = [C+(0)£(0) - B+(0)C(0)]/(2i), A, = C*(0)B(0) - S*(0)C(0). 



The relation of the matrices B(t) and C(t) to the matrix Q{t) and, in view of fl5.q ), to the function 4>t{t) 
yields 
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Statement A. 3 Let the n x n-matrices B(t) and C(t) be solutions to equations in variations (5.4). Then, if 
det C(t) ^ 0, t G [0,T], the matrix Q(t) = i3(i)C _1 (i) satisfies the Riccati matrix equation (5.9). 



Actually, in view of 

and since from C- x (£)C(t) = I follows C , - 1 (t)C(t) + C , - 1 (t)(7(t) = 0, we have 

Q = S(i)C- 1 (t) + B(i)C'- 1 (t) = B(t)(7- 1 (i)-Q(t)C(i)C'- 1 (t) = 

= [-flx„(*)B - ^M^cr 1 - Q[8 w (t).B + vwqc- 1 = 

= -$)xp(t)Q - fira(t) - Qfi pp (t)Q - QSj px (t), 

Q.E.D. 

A similar property is also valid for the matrix 

- Q- 1 + Q-^xxftJQ -1 + tf,«(t)Q~ x + Q -1 iM*) + fi PP (t) = 0. (A.7) 

Statement A. 4 If at the time zero the matrix Q(t) is symmetrical [Q(0) = Q'(0) at t = 0], it is symmetrical 
at any time t € [0,T] [i.e., = Here, A* denotes the transpose to the matrix A. 



Actually, from (5.9) follows 

Q* + rf xx (t) + fi^WQ* + Q%,(*) + Q'^m 1 = 0, 

since we have 



Hence, the matrix Q (£) satisfies equation (5.9) with the same initial conditions as the matrix Q(t) since, as 
agreed, the matrix Q(0) is symmetrical. The validity of the statement follows from the uniqueness of the solution 
to the Cauchy problem. 

Statement A. 5 The imaginary parts of the matrices Q(i) and can be presented in the form 

Im Q(t) = (C+{t))- l D C-\t), (A.8) 
Im Q-\t) = -(B-^ty+DoB-^t). (A.9) 



Here, the matrix Do is defined by relationship (A. 5). 
Actually, by definition we have 

Im Q(t) = l -[Q+(t) Q(t)] = \{[B{t)C- l {t)]+ B{t)C~\t)} = 

= l -[C+{t)]^[B+{t)C(t) - C+{t)B{t)]C-\t) = [C + (t)]- l D C-\t), 



Q.E.D. Similarly, relationship (A.9) can be proved 



Statement A. 6 Let the matrix Q(t) be definite and symmetrical and the components of the vector yUt), 



j = l,7i be the row elements of the matrix C 1 (t) (A. 24). Then the vectors yUt) satisfy the set of equations 



y = -H xp (t)y - Q(t)H pp (t)y. (A.10) 

Actually, we have 
and hence 

C'- 1 = -C- 1 [H vp {t)Q{t) + Hp X {t)]. (A.ll) 



Transposing relationship (All) for the vectors y(t) ( A. 24 ), we obtain equation ( A.1C ), Q.E.D. 
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Remark A. 2 If the matrix 

Q- 1 (t) = C(t)B- 1 (t) 
is definite and symmetrical, then the matrix B^ 1 satisfies the equation 

B- 1 = B- l [H xx {t)Q- x {t)+H xp {t)}. 



The proof is similar to that of Statement A. 6 



(A.12) 



Statement A. 7 If the matrix Dq ( |A.5| ) is positive definite, the relation 

2iB-\t)H xx (t)(B-\t)y = ^[D^B+^it))- 1 } 

is valid. 

Actually, from ( |A.9| ) follows 

S _1 (t)H xx (t) (S _1 (t))* = 

= l -D- l B+{t)[Q-\t) - {Q*{t))- x ]H xx {t){B-\t)Y = 

= l -D^[B-\t){n xx (t)c{t)B-\t)-n xx {t)c + (t){B-\t)Y + 

+H px (t)B(t)B- 1 (t)-H px {t)B*{t)(B- 1 (t))*)(2iIm Q-\t))B* (t)!^ 1 ]' 
= ii?o 1 [S- 1 WB^B-WB* (t) -B-\t)B* (f )] * = 



2^ 



Q.E.D. 



Statement A. 8 If for equations in variations (A..2) the Cauchy problem is formulated as 

B(t)\ t=s = B , C(t)\ t=s =0, Bl=B , 

then the relation 



(A.13) 



(A.14) 



B~ 1 (T)H xx (t)(B~ 1 (T)) t dT = (B - 1 )*A 2 (At,* )A4 1 (Ai,vI/ )(i?- 1 ) t 



(A.15) 



is valid. Here, A&(At, ^o), fc = 1,4 denote the n x n matrices being blocks of the matriciant of the set of 
equations in variations (5.4) 



A(t, * ) 



A*(t,* ) A|(t,* ) 
A|(t,tf ) A^o) 



A(0,* ) =l2nx2n- 



Let us consider an auxiliary Cauchy problem formulated as 

B(t,e)\ t = s = B , C(t,e)\ t = s = el, I=\\S k J nxn . 

Obviously, we have 



(A.16) 



(A.17) 



and 



]im B(t,e) = B(t), lim C(t, e) = C(t) 



Z>o(e) = jtCBo -S *). 
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We assume that the matrix -Do(e) is symmetrical and positive definite for e 7^ 0. Hence, we may use relationship 



(A. 15) and then obtain 



J B-\t, e)H xx (T)(B-\T, e)) l dT = 

S 

-kB -B* )- l B+(T,e)(B-HT,e)y\i. (A.18) 



In view of (A. 16), we have 
and, hence, 



Then we obtain 



B(t, e) = X\(At)B — eA^ (At), 
B-%e) = (l + eJ B - 1 (A 4 - 1 (At)) t ^(At)) J B - 1 (A 4 - 1 (Ai))*+O( e 2 ). 



lim V(r, e)(B-\t, e )) ( |< = [ limjB-^t, e)B*(t, e) B^B*]} ' 



= hmi{S - 1 i3*- e i? -\A 4 - 1 (Ai)) t A*(At)(l-i? - 1 i?*)-i? - 1 i? *+O( e 2 )]} t = 
e— >o e 

= -(B B*)(B^) t X 2 (At)X^(At)(B^) t . 
Substitution of the obtained expression into ( |A. 18 ) yields ( A.15| ). 



Statement A. 9 If the matrix Dq ( |A.5| ) is positive definite and symmetrical and the matrix D$ (A. 6) is zero, 
the following relationships are valid: 

C*(t)D^B t (t)-C(t)D^B + (t) = 

= B{t)D l C+{t) - B'tyDo 1 *?® = 2iI nXn , (A.19) 

C*(t)D - l C t (t)-C(t)D - 1 C+(t) = 

= B{t)D^B+{t) - B*(t)D^ 1 B t (t) = nxn . (A.20) 

Let us consider an auxiliary matrix T(t) of dimension 2n x 2n 



T(t) J- ( D o 1/2 C\t) -D-^B^t) , 
( } _ V2 [ D-^C+it) -D-^B+it) 1 (A ' 21) 



and find its inverse matrix. Direct checking makes us convinced that 

i_ ( -(B(t)D Q 1/2 Y B(, 

V2 { -(c(t)D- 1/2 y c(t)D- 



Actually, we have 

T(t)T- l (t) = 

1 f-[C\t)B*{t) - B\t)C*{t)](D^Y [C\t)B(t) - B\t)C(t)]D^ 



2\-[C+{t)B*{t) - B+(t)C*(t)](D^)* [C+(t)B(t) B+(t)C(t)]D^ 



= _l -D 1/z (2zD )*(D- 1/ y D 1/Z D D~ 

2 { -D- 1/2 D* (D^ 1/2 y l D- 1/2 D D Q 

From the uniqueness of the inverse matrix follows 

i-1/.A rr-l 



H2rix2n- 



T(t)T-\t) = T-\t)T(t) = I 2nx2 „, 
T(t)T _1 (t) = 
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-(BD^yC 1 + BD^C + (BD^ 1 )*B t + BD^B^ 

= ^2nx2n- 



(A.23) 



However, as assigned, we have Dq = Dq, and from Definition (AJ5) follows 

D+ = -^{C + B - B+C)+ = -^{B + C - C+B) = D . 

We then have = D and, hence, (D^ 1 )* = Dq 1 . Then from ( |A.23[ ) we obtain ( |A.19| ) and ( |A.2C| ). 

The following properties of the solutions to the set of equations in variations are dramatically important for 
the construction of quasi-classical asymptotics in the class of functions V\,{Z(t, K), S(t, h)). 

Lemma A.l Let the matrix Dq be diagonal and positive definite and detC(t) =/= 0. The matrix Im Q(t) will 
then be positive definite as well. 

Proof. Let we have Dq = diag (ai, . . . , a n ), a.j > 0, j = 1, n. Denote by the rows of the matrix C _1 (t): 



c-\t) 



( m \ 

m 



(A.24) 



V &(*) J 



Then for an arbitrary complex vector \p\ ^ we obtain 



(A.25) 



The above inequality is true since |j7j(t)| ^ 0, and ay > 0, j = l,n. From this inequality, in view of the 
arbitrariness of the vector p E C™, \p\ ^ 0, follows the validity of the lemma statement. 

Lemma A.2 The matrixC(t) is nondegenerate [det C(t) ^ 0] if the matrix D = (2i)- 1 (C+(0)B(Q)-B+(Q)C(0)) 
is positive definite. 

Proof. Let, for some ti, detC(ti) = be valid. Then a vector k, \k\ ^ exists, such that 

C(ti)-k = 0, (k + C + (h) = 0). 



Since relationship ( A. 6 ) is valid for any t, then 

k+D k = fc+{i[.B + (ti)<7(fi) - C+(t 1 )B(ti)]}fc = 0. 

As agreed, the matrix D is positive definite, and, hence, the above equality holds only for \k\ =0. The 
obtained contradiction proves the lemma. 

Lemma A. 3 (Liouville's lemma) If the matrix Q(t) is continuous, the relation 

t 



exp { - i J Sp [Sj pp {t)Q{t) + ft px (t)] dt} 



/detC(O) 
det C{t) ■ 



(A.26) 



is valid. 



Proof. From (A.2) follows 

C = [S!>„(t)Q(t) + ft px (t)]C, (A.27) 

where the matrix Q(t) is a solution of equation (|5.9| ) and the matrices Sj pp (t) and S) px (t) are continuous. Then 
for the matrix C(t) the Jacobi identity 

t 

detC(t) = [detC(0)]ex P J Sp [Sj pp (t)Q(t) + S) px (t)} dt 



is valid. Raising the left and right parts of the equality to the power —1/2 yields ( A.2q), Q.E.D. 
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